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INTRODUCTION 


Analytical  techniques  exist  for  determining  the  electromagnetic  (EMP) 
fields  that  are  transmitted  into  the  ground  for  given  parameters  of  the  ground 
and  the  incident  field.  The  main  objective  of  this  paper  is  to  indirectly 
measure  the  conductivity  o,  dielectric  constant  er,  and  electric  (E)  field 
below  ground  due  to  an  incident  EMP  field  as  produced  by  the  repetitive  EMP 
simulator  (REPS).  REPS  is  a  horizontal  dipole  radiator  driven  by  a  1-MV 
repetitive  pulse  generator.  The  measurements  were  taken  at  the  Harry  Diamond 
Laboratories  (HDL),  Woodbridge  Research  Facility  (WRF ) ,  Woodbridge,  VA. 

There  are  two  methods  to  accomplish  this  objective: 

(1)  The  determination  of  the  "calculated"  transmitted  E  field,  Et,  from 
the  measured  magnetic  field,  Hx,  at  1  m  above  ground,  and  the  associated 
Maxwell  equations  and  Fresnel  coefficients  in  a  continuous  air/ground  inter¬ 
face.  The  ground  can  be  treated  as  a  good  conductor  whose  most  important 
electrical  parameters  are  conductivity  and  dielectric  permittivity,  e. 

(2)  The  determination  of  the  "measured"  transmitted  E  field  from  the 
induced  voltage,  V,  across  a  buried  capacitive  parallel-plate  E-field  sensor 
with  a  plate  separation  f.  This  V  is  "unfolded"  from  the  measured  voltage, 
VQ,  at  the  sensor  load  through  the  use  of  the  time-domain  and  frequency-domain 
solution  techniques.  The  time-domain  solution  is  derived  from  an  equivalent 
circuit  model  of  the  E-field  sensor.  From  this  solution,  the  sensor  can  be 
characterized  as  an  E-field  sensor  and  an  E  (first  derivative  of  E)  field 
sensor.  The  frequency-domain  solution  technique  depicts  the  behavior  of  the  E 
field  below  ground  through  parametric  variations  of  frequency-independent 
(constant)  ground  parameters.  The  frequency-domain  solution  of  the  same 
equivalent  circuit  model  of  the  E-field  sensor  uses  ground  parameters  that  are 
either  frequency  independent  (constant)  or  frequency  dependent  to  describe  the 
transfer  function.  A,  or  sensor  calibration  of  the  sensor. 

These  two  methods  independently  arrive  at  the  solution  of  the  transmitted 
E  field  below  ground  but  both  depend  on  o  and  er.  The  parameters  o  and  er 
were  the  only  ones  adjusted  to  provide  agreement  between  the  calculated  and 
measured  E  fields  transmitted  below  ground.  When  the  calculated  and  measured 
E  fields  are  in  good  agreement  for  given  o  and  er,  a  conclusion  can  be  drawn 
from  the  results. 

A  flow  chart  of  an  indirect  measure  of  below-ground  E  field,  conductivity, 
and  dielectric  constant  is  shown  in  figure  1. 

This  report  presents  comparisons  between  calculated  and  measured  trans¬ 
mitted  E  fields  using  both  constant  and  frequency-dependent  ground  parameters. 
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Figure  1 .  Flow  chart  of  an  indirect  measure  of  below-ground 
E  field,  conductivity,  and  dielectric  constant. 

2.  ANALYTICAL  CALCULATION  OF  TRANSMITTED  E  FIELD 

When  an  EMP  is  incident  at  the  plane  boundary  of  a  linear,  homogeneous, 
isotropic,  and  conducting  medium,  some  of  it  is  reflected  and  some  is  trans¬ 
mitted  into  the  medium.  The  transmitted  E  field  can  be  found  with  the  use  of 
Fresnel  reflection  and  transmission  coefficients.1-3  The  Fresnel  coefficients 
are  a  function  of  the  electrical  properties  of  the  ground  and  the  incident 
angle  of  the  electromagnetic  wave.  It  is  assumed  in  this  study  that  the 
incident  wave  is  a  linearly  polarized  plane  wave  (of  constant  amplitude  and 
phase)  and  the  air/ground  boundary  is  a  semi -inf ini te  plane. 

The  pertinent  equations  involve  plane  monochromatic  waves  (i.e.,  with 
single  frequency)  as  directly  derived  from  Maxwell's  equations.  Detailed 
derivations  governing  these  equations  can  be  found  in  works  cited  in  the 
Selected  Bibliography.  Figure  2(a)  diagrams  the  wave  vectors  of  the  incident, 
reflected,  and  transmitted  waves  used  in  this  study.  Figure  2(b)  shows  the 
conventional  directions  of  electric  and  magnetic  fields  for  horizontal  polar¬ 
ization. 

C.  Jordan  and  K.  G.  Balmain,  Electromagnetic  Waves  and  Radiating 
Systems,  Prentice-Hall,  Inc.,  Englewood  Cliffs,  NJ  (1968),  2nd  ed./  ch  5,  p 
144  ff . 

2J.  D.  Jackson,  Classical  Electrodynamics,  John  Wiley  s  Sons,  Inc,,  New 
York  (1962),  ch  7,  p  216  ff. 

3 AT.  Born  and  B.  Wolf,  Principles  of  Optics,  Pergamon  Press,  Oxford  (1970), 
fourth  ed.f  ch  1,  p  40;  ch  13,  p  615  ff. 
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Figure  2.  Waves  and  fields:  (a)  wave  vectors  of  incident  reflected  and 
transmitted  waves  and  (b)  conventional  directions  of  electric  and  magnetic 
fields  for  horizontal  polarization. 

The  solution  for  the  incident  electric  field  above  ground  treated  by  Marx4 
is 

Eo(0>)  =  {z0Hx(u>)/[sin  i|>  (l  -  Rh(u)e“^wtd)J}  | sin(utd)/(a)td)  | 2  .  (1) 

E0(u>)  is  a  function  of  the  free-space  wave  impedance  ZQ,  the  magnetic  field 
Hx(u>),  incident  angle  ty,  Fourier  transform  variable  u>  (u  =  2nf),  the  Fresnel 
coefficient  for  horizontal  polarization  Rh(o),  and  time  delay  td.  The  time 
delay  describes  the  difference  of  arrival  time  between  the  incident  and  re¬ 
flected  pulses  at  the  field  measurement  point  above  ground  as  shown  in  figure 
2(a).  The  last  term  of  equation  (1)  is  a  filter  function  that  removes  the 
singularities  at 


<Utd  a  kW 


1,  3,  5,  ... 


Rh(w)  is  found  to  be 


Rh(<*>)  = 


sin  -  [er  -  j(o/e0u>)  -  cos2!);]1/2 
sin  +  [er  -  j(o/e0w)  -  cos2^]1/2 


where  e  is  the  dielectric  permittivity  of  free  space.  The  time  delay  is 


2h  sin 


where  h  is  the  height  of  the  H-field  sensor  above  ground  (1  m)  and  c  is  the 
speed  of  light. 

**Egon  Marx,  Simulator  Fields  and  Ground  Constants/  Harry  Diamond  Labora¬ 
tories,  HDL-TR-1785  (February  1977). 


(4) 


The  transmitted  electric  field  Et(m)  is 

Et(u))  =  Th(uj)e~YdE0(aj)  , 

where  d  is  the  distance  from  the  interface  in  the  soil  and  y  is  the  propaga¬ 
tion  constant,  defined  as 

Y  =  (-u)2pe  +  jonio)1/2  .  (5) 

The  transmission  coefficient5  Th(u>)  is 


Thfio)  = 


2  sin 


{sin  i|»  +  [er  -  j(o/uje0)  -  cos2^]1/2} 


The  time-domain  electric  field  Et(t)  is  numerically  computed  through  an  in¬ 
verse  Fourier  transform6  of  equation  (4).  Finally,  Et(t)  is  averaged  over  a 
12-in.  depth  from  1  to  13  in.  below  the  surface;  the  result  is  taken  as  the  E 
field  at  7  in.  below  the  ground. 


3.  EXPERIMENTALLY  MEASURED  FIELDS 

Field  measurements  were  made  at  the  REPS  facility  at  the  following  loca¬ 
tions  (see  fig.  3): 

(a)  test  point  4  (TP4) ,  close  to  the  centerline  at  x  >  800  ft  and  y  »  82. 5 
ft  south  of  the  centerline,  and 

(b)  test  point  1  (TP1),  off  the  centerline  at  x  *  800  ft  and  y  *  609  ft 
north  of  the  centerline. 

At  each  test  point,  two  field  measurements  were  taken:  (1)  the  total 

magnetic  field,  Hx(t),  at  1  m  above  ground  and  (2)  the  transmitted  component 
of  the  tangential  E  field  averaged  over  a  12-in.  depth  from  1  to  13  in.  below 
the  surface. 

The  Hx(t)  was  measured  with  a  conventional  Stanford  Research  Institute 
(SRI)  cubical  sensor  box.7  Figure  4  shows  the  measured  Hx(t)  at  TP1  and  TP4. 


5Egon  Marx,  Reflected  and  Transmitted  Fields  for  a  Plane-Have  Pulse  Inci¬ 
dent  on  Conducting  Ground,  Harr y  Diamond  Laboratories,  HDL-TR-1740  ( April 
1976). 

6 Alfred  G.  Brandstein  and  Egon  Marx,  Numerical  Fourier  Transform,  Harry 
Diamond  Laboratories,  HDL-TR-1748  ( September  1976). 

7B.  C.  Tupper,  R.  H.  St e hie,  and  R.  T.  Wolfram,  BMP  Instrumentation  Devel¬ 
opment,  Stanford  Research  Institute,  report  7990,  under  contract  to  Harry 
Diamond  Laboratories ,  Contract  DAAK-02-69-C-0674. 
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Figure  3.  Isometric  view  of  buried  E-field  sensor 
(dashed  lines),  measuring  instruments,  and  REPS. 


The  transmitted  E  field  was  measured  by  the  use  of  two  buried  parallel 
aluminum  plates.  These  plates  are  12  in.  long  x  12  in.  wide  x  1/4  in.  thick. 
They  were  separated  a  distance  i  of  1.75  in.  at  TP4  and  2.5  in.  at  TP1 .  In 
both  cases  the  plates  were  inserted  to  achieve  intimate  ground  contact.  The 
voltage  VQ(t)  generated  across  the  plates  by  the  field  was  measured  by  the  use 
of  a  fiber-optic  transmitter  attached  to  the  plate  by  an  R6214,  50-s)  coaxial 
double-shielded  cable.  This  cable  was  1.75  ft  long  at  TP4  and  3  ft  long  at 
TP1 •  The  50-fi  fiber-optic  transmitter  was  connected  to  the  instrumentation 
van  remote -reading  equipment  by  a  fiber-optic  cable  and  a  50-fl  fiber-optic 
receiver.  The  fiber-optic  transmitter  and  receiver  data  link  were  designed 
and  built  by  Jim  Blackburn  of  HDL.  Figure  5  shows  the  measured  sensor  voltage 
VQ(t)  at  TP1  and  TP4. 


I*: 


4.  EFFECTS  OF  CONSTANT  GROUND  PARAMETERS — SENSOR  CHARACTERIZATION 


The  E-field  data  are  obtained  from  the  measured  voltage  by  the  use  of  the 
equivalent  circuit  model  as  shown  in  figure  6.  This  model  is  described  by 
Baum8  as  a  short  dipole  antenna  model.  When  time  variations  are  slow  enough 
that  the  short  antenna  approximation  is  valid,  and  assuming  that  the  edge 
effects  of  the  plates  are  not  a  significant  factor,  the  equivalent  circuit  is 
used  to  represent  the  following  relationship  between  the  voltage  through  the 
load,  VQ(t),  and  the  magnitude  of  the  electric  field,  Et(t): 

/  d  \r  i  VQ(t)  dVQ(t) 

(cs  aT  *  G)(v(t>  ~  vo<*>]  *  — 5^-  ♦  Cc  -at"  '  for  t  >  0  .  (7) 

Here  V(t)  =  Et(tU,  G  *  1/Rg  =  Cga/e,  R^  is  the  load  resistance,  Cg  =  area, 
e/Z  is  the  sensor  capacitance,  and  Cc  is  the  cable  capacitance.  This  is  the 
same  model  for  the  E-field  sensors  used  in  Aurora  with  time-varying  air  con¬ 
ductivity.  9 


v- 


Figure  6.  Equivalent  circuit  of  E- 
field  sensor. 


The  solution  of  differential  equation  (7)  is  composed  of  a  complementary 
and  a  particular  integral.  In  network  terminology,  these  are  also  referred  to 
as  the  natural,  source-free,  or  transient  response  and  the  forced  or  steady- 
state  response,  respectively.  Examination  of  equation  (7)  reveals  that  when 
Rl  is  large  and  Cc  is  zero,  the  source  terms  are  zero  and  the  solution  is  only 
complementary.  Solving  equation  (7)  by  the  method  of  variation  of  parameters, 
the  solution  for  Efc(t)  is 


8C.  E.  Baum,  Electromagnetic  Pulse  Sensor  and  Simulation  Notes,  Vol.  1,  Air 
Force  Weapons  Laboratory,  Note  13  (June  1970), 

9 Rolando  P.  Manriquez,  George  Merkel,  William  D,  Scharf ,  and  Daniel  Spohn, 
Electrically  short  Monopole  Antenna  Response  in  an  Ionized  Air  Environment, 
Determination  of  Ionized  Air  Conductivity,  IEEE  Trans,  Nucl,  Sci,  NS-26 ,  6 
(December  1979),  5012-5018, 


Et(t) 


V«(t)  + 


Cg  -  GR^Cc  -t/RgCg 

J  0 

arlc| 


-t’/RgCo 

x  /  e  VQ(t*)  dt'  ,  for  t  >  0  .  (8) 

0 

Several  important  observations  may  be  made  from  equation  (8).  When  is 
large  and  Cc  is  zero,  and  the  transit  time  (£/c)  of  the  antenna  model  is  long 
compared  to  the  rise  time  x  of  the  incident  pulse  (i.e.,  R^Cg  >>  x),  the  first 
term  on  the  right-hand  side  of  equation  (8)  dominates.  In  this  case,  the 
sensor  can  be  regarded  as  an  E-field  sensor.  On  the  other  hand,  when  R^Cg  << 
x,  the  second  term  of  equation  (8)  dominates  and  the  sensor  can  be  regarded  as 
an  E-fielda sensor.  Otherwise,  the  sensor  can  be  described  as  a  combination  of 
an  E-  and  E-field  sensor.  The  transmitted  electric  field  Et(t)  can  be  numeri¬ 
cally  computed  from  equation  (8)  or,  alternatively,  the  voltage  V(t)  *  £Et(t) 
can  be  obtained  by  solving  the  differential  equation  (7)  by  a  Runge-Kutta  or 
Gear  method.  The  solutions  of  the  first  and  second  terms  of  equation  (8)  are 
shown  in  figure  7  at  TP1  and  figure  8  at  TP4,  respectively,  with  a  constant  a 
=  0.007  mho/m  and  e_  =  15. 


THE  (I)  (xIO  *) 

Figure  7.  Results  of  first  term  (solid  line)  and  second  term 
(dashed  line)  of  equation  (8)  at  TP1  with  o  -  0.007  aho/m  and 
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Figure  8.  Results  of  first  term  (solid  line)  and  second  term 
(dashed  line)  of  equation  (8)  at  TP4  with  o  «  0.007  mho/m  and 


Actual  measurements  of  a  and  er  were  not  available  for  the  time  this  test 
was  performed.  However*  previous  data  collected  by  the  National  Bureau  of 
Standards  (NBS)  show  that  the  ground  conductivity  is  approximately  0.007  mho/m 
and  the  er  is  15  at  1  MHz.  These  data  are  discussed  elsewhere10  and  measured 
for  a  limited  frequency  range.  The  results  from  equation  (8)  and  the  time- 
domain  Fourier  transform  of  equation  (4)  are  shown  in  figures  9  to  13,  at  TP4, 
for  o  varied  with  0.001,  0.007,  and  0.02  mho/m  at  er  =  15,  and  er  varied  with 
1,  15,  and  80  at  0  >  0.007  mho/m.  Figure  14  shows  the  comparison  between  the 
results  of  equation  (8)  and  the  time-domain  Fourier  transform  of  equation  (4), 
at  TP1,  for  a  -  0.007  mho/m  and  er  *  15.  The  significance  of  the  parametric 
variational  effects  to  the  expected  values  at  the  extreme  is  apparent.  As  o 
and  er  increase,  the  amplitude  of  the  electric  field  decreases.  The 
waveshapes  at  late  times  and  low  frequencies  are  somewhat  altered  at  higher 
conductivities.  The  peak  amplitude  is  particularly  sensitive  to  the  changes  of 
the  dielectric  constant  at  higher  frequencies. 


10JToroan  V.  Hill,  Effect  of  Frequency-Dependent  Soil  Parameters  on  Reflec¬ 
tion  Coefficients,  Harry  Diamond  Laboratories,  HDL-TR-2004  (December  1982), 


17 


TIME  (t)  (xIO  ) 

Figure  13.  Comparison  between  calculated  (dashed  line)  and  measured  (solid 
line)  transmitted  electric  fields  at  TP4,  with  o  *  0.007  mho/m  and  e_  =  80 


One  of  the  unique  features  of  the  measuring  system,  depicted  in  figure  3, 
is  the  fiber-optic  system.  The  advantage  of  implementing  this  system  is  to 
electrically  isolate  the  E-field  sensor  from  the  instrumentation  van,  thereby 
eliminating  the  need  for  a  long  cable  between  the  sensor  and  the  van. 

5.  EFFECTS  OF  FREQUENCY-DEPENDENT  GROUND  PARAMETERS— SENSOR  CALIBRATION 

The  use  of  constant  (frequency-independent)  values  for  o  and  results  in 
a  sensor  calibration  (or  transfer  function)  A  that  is  a  constant.  This  A  can 
be  used  to  determine  the  transmitted  E  field  across  the  parallel-plate  sensor 


Et(w) 


Vq(co) 


However,  in  reality,  o  and  er  are  frequency  dependent,  and  for  larger  varia¬ 
tions  of  frequency,  a  more  accurate  calibration  of  the  buried  E-field  sensor 
must  include  a  frequency-dependent  transfer  function. 

Let  A(tu)  be  the  transfer  function  of  the  buried  E-field  sensor  as  deter¬ 
mined  by  taking  the  Laplace  transform  of  equation  (7)  in  the  s-domain  (s  = 
joi).  In  general,  the  transfer  function  is  a  complex  quantity  and  can  be 
written  as 

A(u>)  -  V(0))/Vo(aj)  .  (10) 

A(oj)  is  also  stated  in  terms  of  magnitude  and  phase  as 


A(<u )  =  |A(a,)|ej*(u,)  , 


where  | A (u> )  |  is  the  amplitude-response  function  and  *(u>)  is  the  phase-shift 
function  of  the  sensor.  The  transfer  function  depends  on  the  circuit  parame¬ 
ters  as 

M  +  SB 

A(“} -rr^r  '  (12> 


where  M  *  (l  +  R^G ) /R^  and  B  »  Cg  +  Cc.  The  amplitude-response  function 
|A(w)|  is 

|A(u)|  -  [l/(G2  +  <02c|)][(mG  +  u>2BC8)2  +  0)2(bG  -  CgM)2]1/2  ,  (13) 

and  the  phase-shift  function  $(b>)  is 


♦(w)  «  tan 


oi(BG  -  CBM) 
MG  +  U)2BC8 


The  results  for  |Vo((0)|,  |v(uj)|,  |A(u>)|,  and  4(u>)  as  a  function  of  fre¬ 

quency  are  shown  in  figures  IS  and  16  at  TP4,  respectively,  for  a  constant  o  ■ 
0.007  mho/m  and  ef  *  IS.  Ideally,  this  sensor  should  produce  an  amplitude 
frequency  response  that  looks  "flat”  in  the  frequency  band  of  interest  and  a 
phase  shift  that  is  a  linear  function  of  frequency.  In  other  words,  the 
spectrum  of  the  measured  input  voltage  V(ui)  is  identical  to  the  spectrum  of 
the  output  voltage  vo(<o)  as  expressed  in  equation  (10).  This  means  that  the 
input  voltage  is  passed  undistorted  by  the  measuring  system.  But  for  some 
cases,  when  the  amplitude  and  phase  frequency  response  are  functions  of  a  and 
er  that  vary  with  frequency  and  moisture  content,  the  output  voltage  may  be 
substantially  different  from  the  input.  From  these  viewpoints,  depending  on 
the  ground  parameters,  the  sensor's  transfer  function  could  appreciably  alter 
or  distort  the  output  voltage. 

Finally,  the  "unfolded"  measured  transmitted  electric  field  is 

Et(oj)  =  A(u))Vq(o >)/Z  .  (15) 


FREQUENCY  (Hi) 

Figure  15.  Magnitude  of  transfer  function  A(u)  (solid  line), 
measured  sensor  voltage  VQ(u)  (dashed  line),  and  "unfolded" 
sensor  voltage  V(u)  (dash-dot  line)  at  TP4,  with  o  *  0.007 
nho/m  and  er  *  15. 
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Figure  16.  Phase  shift  function  $(u>)  response  of  sensor  at  TP4, 
with  o  *>  0.007  mho/m  and  er  =  15. 


Several  studies  have  been  conducted  concerning  the  measurements  and  theo¬ 
retical  formulations  of  the  electrical  properties  of  the  soil,  namely,  a  and 
er  as  functions  of  frequency  and  moisture  content.  Longmire  and  Smith11 
developed  a  universal  formula  for  a  and  er  over  the  frequency  range  of  5  Hz  to 
3  *  1012  -Hz,  based  on  Scott's  data12  for  soils  and  WilXenfeld's  data  for  some 
concrete  and  grout  samples  (Wilkenf eld's  data  can  be  found  in  Longmire  and 
Smith1 1 ) . 


nC.  L,  Longmire  end  K.  S,  Smith,  A  Universal  Impedance  for  Soils,  Mission 
Research  Corp.,  Santa  Barbara,  CA,  Contract  Mo,  DNA00 1-7 5-C-0094  (October 
1975), 

12j.  H.  Scott,  electrical  and  Magnetic  Properties  of  Rock  and  Soils,  Mote 
18,  electromagnetic  Pulse  Theoretical  Motes,  Air  rorce  Weapons  Laboratory,  BMP 
2-1  (April  1971) i  also  U,S.  Geological  Survey  Technical  Letter,  Special  Proj¬ 
ect  16  (26  May  1966), 
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The  0  and  er  derived  from  Longmire  and  Smith's  "universal  RC  network 
model"  are 


:r  *  e»  +  I  - T2 - T7 


(relative)  , 


0*0.  +  2  ire.  V  a_f„ 
o  o  L  n  n 


<-*, . 


where 


N  =  13, 
e«»  *  5* 

f  *  frequency  (Hz), 

an  *  the  constant  coefficients  (see  table  1), 
fn  *  F(P)fn(10%), 
fn(1°*)  *  HZf 

F(P)  *  (P/10)1,28, 

P  *  water  content  (percent),  and 
oQ  *  8.0  x  io“3  (P/10)1,54  (mho/m). 

Figure  17  shows  the  ground  conductivity  a  versus  frequency  for  various  volume 
percentages  of  water.  Figure  18  shows  the  dielectric  constant  er  versus 
frequency  for  various  volume  percentages  of  water. 


TABLE  1.  COEFFICIENT  *n  FOR  UNIVERSAL  SOIL 
(see  eq  (16)  and  (17)) 
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an 

n 

an 
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6 
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11 

9.80 
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7 
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12 

3.92 
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3 
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8 
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13 

1.73 
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4 

3.38  x  103 

9 
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5 

5.26  x  in* 

10 

2.17 
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Because  of  the  unavailability  of  Woodbridge's  soil  data  over  a  wide  range 
of  frequency,  it  was  necessary  to  implement  Longmire  and  Smith's  universal 
formula  for  o  and  er  in  the  program.  However,  some  old  data  taken  by  NBS  for 
Woodbridge's  soil  show  relatively  low  a  and  er;  Hill10  discusses  these  data. 
In  the  analysis,  10-percent  soil  moisture  content,  a  *  0.007  mho/m,  and  er  * 
15  provided  close  agreement  between  calculated  and  measured  transmitted  E 
fields  below  ground.  The  results  for  j  v0 < <*> >  j  *  |v(w)|,  |a(u>)|,  and  with 
10-percent  moisture  content  at  TP4  are  shown  in  figures  19  and  20.  The  com¬ 
parison  between  the  inverse  Pourier  transform  of  the  calculated  transmitted 
electric  field  (eq  (4))  and  the  "unfolded"  measured  transmitted  electric  field 
(eq  (15))  using  a  and  er  dependent  with  frequency  at  10-  and  25-percent  mois¬ 
ture  content  are  shown  in  figures  21  and  22  at  TP4,  and  figures  23  and  24  at 
TP1.  Indeed,  the  transfer  function  of  the  sensor  is  highly  sensitive  to  the 
electrical  parameters  of  the  soil. 


A-5V.,  B-10%,  C-15%,  D-20'/«,  E-25% 
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Figure  17.  Conductivity  versus  frequency  for  various  volume 
percentages  of  water. 


10  Norman  V .  Mill,  Effect  of  Frequency-Dependent  Soil  Parameters  on  Reflec¬ 
tion  Coe f fide.,  j,  Harry  Diamond  laboratories,  HDL-TR-2004  (December  1982), 


FREQUENCY  (Hz) 

Figure  20.  Phase  shift  function  $(o>)  response  of  sensor  at  TP4,  using 
Longmire's  soil  data  ( 1 0-percent  moisture  content). 


Figure  21.  Comparison  between  calculated  (dashed  line)  and  measured  (solid 
line)  transmitted  electric  fields  at  TP4,  using  Longmire's  soil  data  (10- 
percent  moisture  content). 
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Figure  24.  Comparison  between  calculated  (dashed  line)  and 
measured  (solid  line)  transmitted  electric  fields  at  TP1 ,  using 
Longmire's  soil  data  (25-percent  moisture  content). 

6.  CONCLUSION  AND  RECOMMENDATIONS 

This  paper  documents  an  attempt  to  experimentally  measure  the  E-field 
component  of  an  EMP  below  ground  and  to  compare  the  results  to  an  analytical 
calculation.  The  results  of  equations  (4)  and  (15)  showed  good  agreement  for 
10-percent  moisture  content,  constant  o  =  0.007  mho/m,  and  er  *  15. 

The  measurements  can  be  improved  by  the  use  of  a  differential-mode  voltage 
probe  to  measure  the  transmitted  E  field.  The  RG214  cable  can  be  removed  and 
the  sensor  directly  connected  to  the  fiber-optic  transmitter.  Another  im¬ 
provement  would  be  to  accurately  determine  o  and  e  over  a  wide  range  of 
frequency  and  depth  at  the  same  location  where  the  fields  were  measured,  and 
at  about  the  same  time.  The  availability  of  more  soil  data  would  reduce  the 
uncertainty  in  the  sensor  calibration. 

The  assumption  that  the  EMP  was  a  plane  wave  over  a  homogeneous  plane 
semi -infinite  ground  in  the  far-field  radiation  zone  may  be  justified  by  the 
quality  of  the  results.  The  transmitted  E  field  vanishes  at  late  times  (>1 
lis)  but  the  transmitted  H  fields  may  not.  Also  not  taken  into  account  were 
the  multiple  reflection  of  the  fields  and  the  effects  of  dispersion  due  to  the 
existence  of  different  layers  of  strata  below  ground.  These  uncertainties  can 
be  resolved  by  measurement  of  the  E  and  H  fields  at  different  depths  below  the 
air/ground  interface. 


I 


Analysis  has  shown  that  it  may  be  possible  to  indirectly  measure  o  and  er 
as  a  function  of  moisture  content  and  frequency  with  the  parallel-plate  E- 
field  sensor.  The  sensitivity  of  the  sensor  is  demonstrated  through  the 
sensitivity  (S)  analysis  of  the  equivalent  circuit  model  of  the  sensor's 
response  to  varying  moisture  content  (P),  i.e.,  S  =  3V/3P  (see  fig.  21  to  24). 

Perhaps  the  most  significant  aspect  of  this  effort  is  that,  through  an 
adequate  calibration  of  the  buried  E-field  sensor,  a  method  now  exists  for  the 
immediate  and  relatively  easy  indirect  determination  of  ground  conductivity 
and  dielectric  constant.  This  method,  through  the  Fourier  transform,  could 
then  be  made  available  as  frequency-domain  data  and  applied  as  derived  to  all 
EMP  coupling  programs. 

The  curve-fitting  equations  (eq  (16)  and  (17))  used  by  Longmire  and 
Smith11  that  determined  the  er  and  o  based  on  Scott  and  Wilkinfeld's  data  can 
be  further  modified  by  adjusting  the  necessary  coefficients  in  the  equations 
to  obtain  a  closer  correlation  between  the  calculated  and  measured  transmitted 
electric  fields.  A  computer-aided  optimization  procedure13  is  needed  to 
accomplish  this  task.  This  curve-fitting  method  can  analytically  improve  the 
determination  of  the  o  and  £r  for  Woodbridge's  soil. 

The  use  of  calibrated,  shallow,  buried  parallel  plates  should  be  made  part 
of  all  field-test  system  programs  because  it  is  a  simple,  inexpensive  method 
of  determining  the  soil  conductivity  at  the  same  time  that  the  experimental 
coupling  data  are  collected  on  the  system.  Thus,  a  conductivity  measurement 
made  at  the  beginning  of  each  test  day  can  be  used  to  predict  the  signal 
levels  expected.  And,  in  addition,  an  accurate  evaluation  of  experimentally 
collected  data  can  then  be  used  by  the  analyst  to  predict  the  levels  of  in¬ 
duced  signals  for  any  conditions  of  soil. 

Future  efforts  will  be  to  explore  ways  of  improving  the  measurement  method 
by  (1)  the  determination  of  the  effects  of  RG214  cable  on  the  measurement, 

(2)  an  independent  direct  measurement  of  o,  er,  and  P  by  whatever  means, 

(3)  the  use  of  a  different  buried  sensor  (dipole,  magnetic  loop,  two  parallel 
cylinders,  two  parallel  spheres,  etc),  and  (4)  the  fteasurement  of  conduction 
current  density  (oEt)  below  the  ground  using  the  parallel-plate  sensor.  This 
last  measurement  can  be  performed  by  connecting  a  large  resistor  between  the 
sensor  and  a  short  RG214  cable.  In  series  with  the  cable  will  be  an 
impedance-matching  device.  This  device  will  match  the  high-impedance  sensor 
system  (sensor  and  resistor)  to  the  low-impedance  data-link  system  (cable  and 
fiber-optic  system).  The  use  of  an  impedance-matching  device  makes  it  possi¬ 
ble  to  directly  measure  the  induced  sensor  voltage.  The  signal  propagated 
below  the  ground  will  be  produced  by  the  REPS. 


“C.  L.  Longmire  and  K.  S.  Smith,  A  Universal  Impedance  for  Soils,  Mission 
Research  Corp,,  Santa  Barbara,  CA,  Contract  No,  DNA001-75-C-0094  (October 
197 5). 

13*.  Fletcher  and  M,  J,  D,  Powell ,  A  Rapidly  Convergent  Descent  Method  for 
Minimization,  Computer  J,  6  (1963),  163 . 
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